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Abstract 

A binned Dalitz plot analysis of B^ — > DK^ decays, with D — > Ag 7r + 7r~ and D — > 
K®K + K~, is performed to measure the CP- violating observables x± and y± which 
are sensitive to the CKM angle 7. The analysis exploits 1.0 fb _1 of data collected 
by the LHCb experiment. The study makes no model-based assumption on the 
variation of the strong phase of the D decay amplitude over the Dalitz plot, but uses 
measurements of this quantity from CLEO-c as input. The values of the parameters 
are found to be x_ = (0.0±4.3± 1.5±0.6) x 1(T 2 , y_ = (2.7±5.2±0.8±2.3) x 10~ 2 , 
x+ = (-10.3±4.5±1.8±1.4) x 1(T 2 and y + = (-0.9±3.7±0.8±3.0) x 1(T 2 . The 
first, second, and third uncertainties are the statistical, the experimental systematic, 
and the error associated with the precision of the strong-phase parameters measured 
at CLEO-c, respectively. These results correspond to 7 = (44 l 3 g)°, with a second 
solution at 7 — > 7 + 180°, and re = 0.07 ± 0.04, where re is the ratio between the 
suppressed and favoured B decay amplitudes. 
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1 Introduction 



A precise determination of the Unitarity Triangle angle 7 (also denoted as ^3), is an im- 
portant goal in flavour physics. Measurements of this weak phase in tree-level processes 
involving the interference between b — > cus and b — > ucs transitions are expected to 
be insensitive to new physics contributions, thereby providing a Standard Model bench- 
mark against which other observables, more likely to be affected by new physics, can be 
compared. A powerful approach for measuring 7 is to study CP-violating observables 
in B ± — > DK ± decays, where D designates a neutral D meson reconstructed in a final 
state common to both D° and D° decays. Examples of such final states include two-body 
modes, where LHCb has already presented results [j], and self CP-conjugate three-body 
decays, such as K®n + ii~ and K®K + K~ , designated collectively as K®h + h~. 

The proposal to measure 7 with B ± — » DK ± , D — > K®h + h~ decays was first made 
in Refs. [2j[3]. The strategy relies on comparing the distribution of events in the D — > 
K®h + h~ Dalitz plot for B + — > DK + and B~ — > DK~ decays. However, in order to 
determine 7 it is necessary to know how the strong phase of the D decay varies over the 
Dalitz plot. One approach for solving this problem, adopted by BaBar |2]-[6] and Belle [7|- 
M, is to use an amplitude model fitted on flavour-tagged D — > K®h + h~ decays to provide 



this input. An attractive alternative (2,10,11 is to make use of direct measurements 
of the strong phase behaviour in bins of the Dalitz plot, which can be obtained from 
quantum-correlated DD pairs from ^(3770) decays and that are available from CLEO- 



12 , thereby avoiding the need to assign any model-related systematic uncertainty. A 



first model-independent analysis was recently presented by Belle 13 using B — >■ DK ± , 
D —7- K®7t + tt~ decays. In this Letter, pp collision data at = 7 TeV, corresponding 
to an integrated luminosity of 1.0 fb _1 and accumulated by LHCb in 2011, are exploited 
to perform a similar model-independent study of the decay mode B ± — > DK ± with 
D — > A'g7r + 7r~ and D — > K^K + K~ . The results are used to set constraints on the value 
of 7. 



2 Formalism and external inputs 

The amplitude of the decay B + — > DK + , D — > K®h + h~ can be written as the superpo- 
sition of the B + — > D°K + and B + — > D°K + contributions as 

A B (m 2 + , m 2 _)=A + r B e i{5B+ ^A. (1) 

Here m 2 + and m?_ are the invariant masses squared of the K®h + and K®h~ combinations, 
respectively, that define the position of the decay in the Dalitz plot, A = A(m'^_,m 2 _) 
is the D° — > K®h + h~ amplitude, and A = A(m+,m 2 _) the D° — > K®h + h~ amplitude. 
The parameter r#, the ratio of the magnitudes of the B + — > D°K + and B + — > D°K + 



amplitudes, is ~0.1 14 , and S B is the strong-phase difference between them. The equiv- 
alent expression for the charge-conjugated decay B~ — > DK~ is obtained by making the 
substitutions 7 — > —7 and A «-» A. Neglecting CP violation, which is known to be small 



1 



in D° — D° mixing and Cabibbo-favoured D meson decays [15], the conjugate amplitudes 
are related by A(m 2 h ,m 2 _) = A{m 2 _,m 2 + ). 

Following the formalism set out in Ref. [2], the Dalitz plot is partitioned into 2N 
regions symmetric under the exchange m+ -H- m 2 _. The bins are labelled from —N to +iV 
(excluding zero), where the positive bins satisfy m 2 _ > m 2 + . At each point in the Dalitz 
plot, there is a strong-phase difference 8 B (m\,m 2 _) = arg A — argA between the D° and 
D° decay. The cosine of the strong-phase difference averaged in each bin and weighted 
by the absolute decay rate is termed Cj and is given by 

J Vi (\A\\A\ cos 5 D )dV 

(2) 



j Vi \A\ 2 dV ^1 j Vi \A\ 2 dV 

where the integrals are evaluated over the area T> of bin i. An analogous expression may 
be written for s$, which is the sine of the strong-phase difference within bin i, weighted by 
the decay rate. The values of Cj and s, can be determined by assuming a functional form 
for \A\, \A\ and 8d, which may be obtained from an amplitude model fitted to flavour- 
tagged D° decays. Alternatively direct measurements of q and Sj can be used. Such 
measurements have been performed at CLEO-c, exploiting quantum-correlated DD pairs 
produced at the ^>(3770) resonance. This has been done with a double-tagged method in 
which one D meson is reconstructed in a decay to either K®h + h~ or K®h + h~ , and the 
other D meson is reconstructed either in a CP eigenstate or in a decay to K®h + h~ . The 
efficiency-corrected event yields, combined with flavour-tag information, allow q and Si 
to be determined [2 10 11 . The latter approach is attractive as it avoids any assumption 



about the nature of the intermediate resonances which contribute to the K^h + h~ final 
state; such an assumption leads to a systematic uncertainty associated with the variation 
in 5d that is difficult to quantify. Instead, an uncertainty is assigned that is related to 
the precision of the q and Sj measurements. 

The population of each positive (negative) bin in the Dalitz plot arising from B + 
decays is N+ t (iVj^), and that from B~ decays is (ACJ. From Eq. ([!]) it follows that 

N+ = h B+ 
N£ = h B - 



K Ti + (x 2 + + y 2 + )K ±i + 2yjKiK-i{x + c±i =F y+s±i 



K ±i + {xl + yl)K Tl + 2^K l K_ l {x.c ±i ± y_s ±i ) , (3) 



where h B ± are normalisation factors which can, in principle, be different for B + and B~ 
due to the production asymmetries, and Ki is the number of events in bin i of the decay 
of a flavour tagged D° — > K®h + h~ Dalitz plot. The sensitivity to 7 enters through the 
Cartesian parameters 

x± = r B cos(5b ± 7) and y± = r B sm(8 B ± 7). (4) 

In this analysis the observed distribution of candidates over the D — > K®h + h~ Dalitz 
plot is used to fit x±, y± and h B ±. The parameters c$ and Sj are taken from measurements 



performed by CLEO-c 12 . In this manner the analysis avoids any dependence on an 
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amplitude model to describe the variation of the strong phase over the Dalitz plot. A 
model is used, however, to provide the input values for K^. For the D° — y K®tt + tt~ decay 
the model is taken from Ref. (5] and for the D° — y K®K + K~ decay the model is taken 
from Ref. [6] . This choice incurs no significant systematic uncertainty as the models have 
been shown to describe well the intensity distribution of flavour-tagged D° decay data. 

The effect of D°—D° mixing is ignored in the above discussion, and was neglected in the 
CLEO-c measurements of q and Sj as well as in the construction of the amplitude model 
used to calculate Ki. This leads to a bias of the order of 0.2° in the 7 determination 16 
which is negligible for the current analysis. 

The CLEO-c study segments the .?C?7r + 7r~ Dalitz plot into 2x8 bins. Several bin 
definitions are available. Here the 'optimal binning' variant is adopted. In this scheme 
the bins have been chosen to optimise the statistical sensitivity to 7 in the presence of 
a low level of background, which is appropriate for this analysis. The optimisation has 
been performed assuming a strong-phase difference distribution as predicted by the BaBar 
model presented in Ref. [5]. The use of a specific model in defining the bin boundaries 
does not bias the q and measurements. If the model is a poor description of the 
underlying decay the only consequence will be to reduce the statistical sensitivity of the 
7 measurement. 

For the K®K + K~ final state q and S{ measurements are available for the Dalitz plot 
partitioned into 2 x 2, 2 x 3 and 2x4 bins, with the guiding model being that from the 
BaBar study described in Ref. [6]. The bin boundaries divide the Dalitz plot into bins of 
equal size with respect to the strong-phase difference between the D° and D° amplitudes. 
The current analysis adopts the 2x2 option, a decision driven by the size of the signal 
sample. The binning choices for the two decay modes are shown in Fig. [TJ 
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Figure 1: Binning choices for (a) D -y K^-k + tt and (b) D — y K®K + K . The diagonal line 
separates the positive and negative bins. 
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3 The LHCb detector 



The LHCb detector 17 is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5. The detector includes a high precision tracking system consisting of a 
silicon-strip vertex detector surrounding the pp interaction region, a large-area silicon-strip 
detector (VELO) located upstream of a dipole magnet with a bending power of about 
4Tm, and three stations of silicon-strip detectors and straw drift-tubes placed down- 
stream. The combined tracking system has a momentum resolution of (0.4 - 0.6)% in the 
range of 5 - 100 GeV/ c, and an impact parameter (IP) resolution of 20 \xm for tracks with 
high transverse momentum (px)- The dipole magnet can be operated in either polarity 
and this feature is used to reduce systematic effects due to detector asymmetries. In the 
data set considered in this analysis, 58% of data were taken with one polarity and 42% 
with the other. Charged hadrons are identified using two ring-imaging Cherenkov (RICH) 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter 
and a hadronic calorimeter. Muons are identified by a system composed of alternating 
layers of iron and multiwire proportional chambers. 

A two-stage trigger is employed. First a hardware-based decision is taken at a fre- 
quency up to 40 MHz. It accepts high transverse energy clusters in either the electro- 
magnetic calorimeter or hadron calorimeter, or a muon of high px- F° r this analysis, it 
is required that one of the charged final-state tracks forming the B candidate points 
at a deposit in the hadron calorimeter, or that the hardware-trigger decision was taken 
independently of these tracks. A second trigger level, implemented in software, receives 



1 MHz of events and retains ~0.3% of them 18 . It searches for a track with large px and 



large IP with respect to any pp interaction point which is called a primary vertex (PV). 
This track is then required to be part of a two-, three- or four-track secondary vertex with 
a high px sum, significantly displaced from any PV. In order to maximise efficiency at 
an acceptable trigger rate, the displaced vertex is selected with a decision tree algorithm 
that uses px, impact parameter, flight distance and track separation information. Full 
event reconstruction occurs offline, and a loose preselection is applied. 

Approximately three million simulated events for each of the modes B — > 
D(K®n + 'K~)K ± and B — > D(K^tt + tt^)7c ± , and one million simulated events for each of 
B ± ->■ D(K®K + K~)K ± and B ± ->■ D^K+K-)^ are used in the analysis, as well as a 
large inclusive sample of generic B — > DX decays for background studies. These samples 



are generated using a version of Pythia 6.4 19 tuned to model the pp collisions 20 



EvtGen 21 encodes the particle decays in which final state radiation is generated us- 



ing Photos 22 . The interaction of the generated particles with the detector and its 



response are implemented using the Geant4 toolkit 23 as described in Ref. 24 



4 Event selection and invariant mass spectrum fit 

Selection requirements are applied to isolate both B ± — > DK ± and B ± —t Dn^ 1 candi- 
dates, with D — > K®h + h~. Candidates selected in the Cabibbo-favoured B — > 
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decay mode provide an important control sample which is exploited in the analysis. 

A production vertex is assigned to each B candidate. This is the PV for which the 
reconstructed B trajectory has the smallest IP % 2 , where this quantity is defined as the 
difference in the \ 2 fit of the PV with and without the tracks of the considered particle. 
The Kg candidates are formed from two oppositely charged tracks reconstructed in the 
tracking stations, either with associated hits in the VELO detector (long K® candidate) 
or without (downstream K® candidate). The IP \ 2 with respect to the PV of each of 
the long (downstream) K® daughters is required to be greater than 16 (4). The angle 9 
between the candidate momentum and the vector between the decay vertex and the 
PV, expected to be small given the high momentum of the B meson, is required to satisfy 
cosd > 0.99, reducing background from combinations of random tracks. 

The D meson candidates are reconstructed by combining the long (downstream) K® 
candidates with two oppositely charged tracks for which the values of the IP \ 2 with 
respect to the PV are greater than 9 (16). In the case of the D — > K®K + K~ a loose 
particle identification (PID) requirement is placed on the kaons to reduce combinatoric 
backgrounds. The IP x 2 of the candidate D with respect to any PV is demanded to be 
greater than 9 in order to suppress directly produced D mesons, and the angle 9 between 
the D candidate momentum and the vector between the decay and PV is required to 
satisfy the same criterion as for the K® selection (cos^ > 0.99). The invariant mass 
resolution of the signal is 8.7 MeV/c 2 (11.9 MeV/c 2 ) for D mesons reconstructed with 
long (downstream) Kg candidates, and a common window of ±25 MeV/c 2 is imposed 
around the world average D° mass [15]. The K® mass is determined after the addition 
of a constraint that the invariant mass of the two D daughter pions or kaons and the 
two daughter pions have the world average D mass. The invariant mass resolution is 
2.9 MeV/c 2 (4.8 MeV/c 2 ) for long (downstream) K® decays. Candidates are retained for 
which the invariant mass of the two K® daughters lies within ±15 MeV/c 2 of the world 
average K° mass (THj. 

The D meson is combined with a candidate kaon or pion bachelor particle to form the 
B candidate. The IP \ 2 °f t ne bachelor with respect to the PV is required to be greater 
than 25. In order to ensure good discrimination between pions and kaons in the RICH 
system only tracks with momentum less than 100 GeV/c are considered. The bachelor is 
considered as a candidate kaon (pion) according to whether it passes (fails) a cut placed on 
the output of the RICH PID algorithm. The PID information is quantified as a difference 
between the logarithm of the likelihood under the mass hypothesis of a pion or a kaon. 
Criteria are then imposed on the B candidate: that the angle between its momentum and 
the vector between the decay and the PV should have a cosine greater than 0.9999 for 
candidates containing long K® decays (0.99995 for downstream K® decays); that the B 
vertex-separation \ 2 with respect to its PV is greater than 169; and that the B IP x 2 with 
respect to the PV is less than 9. To suppress background from charmless B decays it is 
required that the D vertex lies downstream of the B vertex. In the events with a long Kg 
candidate, a further background arises from B^ 1 — > Dh^, D — > 7T + 7T~h + h~ decays, where 
the two pions are reconstructed as a long K® candidate. This background is removed by 
requiring that the flight significance between the D and K® vertices is greater than 10. 
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In order to obtain the best possible resolution in the Dalitz plot of the D decay, and 
to provide further background suppression, the B, D and K® vertices are refitted with 
additional constraints on the D and K® masses, and the B momentum is required to point 
back to the PV. The x 2 P er degree of freedom of the fit is required to be less than 5. 

Less than 0.4% of the selected events are found to contain two or more candidates. In 
these events only the B candidate with the lowest \ 2 P er degree of freedom from the refit 
is retained for subsequent study. In addition, 0.4% of the candidates are found to have 
been reconstructed such that their D Dalitz plot coordinates lie outside the defined bins, 
and these too are discarded. 

The invariant mass distributions of the selected candidates are shown in Fig. [2] for 
B ± — > DK ± and B ± — > Dir^, with D — > K®TT + 7r~ decays, divided between the long 
and downstream K® categories. Figure [3] shows the corresponding distributions for final 
states with D — > K®K + K~ , here integrated over the two categories. The result of 
an extended, unbinned, maximum likelihood fit to these distributions is superimposed. 
The fit is performed simultaneously for B — > DK ± and B — > Dn^, including both 
D — > K®7t + 7t~ and D — > K®K + K~ decays, allowing several parameters to be different 
for long and downstream K® categories. The fit range is between 5110 MeV/c 2 and 
5800 MeV/c 2 in invariant mass. At this stage in the analysis the fit does not distinguish 
between the different regions of Dalitz plot or B meson charge. The purpose of this 
global fit is to determine the parameters that describe the invariant mass spectrum in 
preparation for the binned fit described in Sect. [5j 

The signal probability density function (PDF) is a Gaussian function with asymmetric 
tails where the unnormalised form is given by 

f exp[-(m - mo) 2 /{2a 2 + a L (m - m ) 2 )],m < m ; , . 
f(m;m ,a L ,a R ,a) = < 5 

^ exp| — [m — mo) / \2o + otR\rn — mo) )\,m > mo; 

where m is the candidate mass, mo the B mass and a, o:l, and an are free parameters in 
the fit. The parameter mo is taken as common for all classes of signal. The parameters 
describing the asymmetric tails are fitted separately for events with long and downstream 
K® categories. The resolution of the Gaussian function is left as a free parameter for the 
two K® categories, but the ratio between this resolution in B — > DK ± and B ± — > 
decays is required to be the same, independent of category. The resolution is determined 
to be around 15 MeV/c 2 for B — > Dn^ decays of both K® classes, and is smaller by a 
factor 0.95 ±0.06 for B — > DK ± . The yield of B — > Dn^ 1 candidates in each category is 
determined in the fit. Instead of fitting the yield of the B — > DK ± candidates separately, 
the ratio 1Z = NlB^ 1 — > DK ± )/N(B ± — > D^) is a free parameter and is common across 
all categories. 

The background has contributions from random track combinations and partially re- 
constructed B decays. The random track combinations are modelled by linear PDFs, 
the parameters of which are floated separately for each class of decay. Partially recon- 
structed backgrounds are described empirically. Studies of simulated events show that 
the partially reconstructed backgrounds are dominated by decays that involve a D meson 
decaying to K®h + h~ . Therefore the same PDF is used to describe these backgrounds as 
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Figure 2: Invariant mass distributions of (a,c) — > DK^ and (b,d) — > Dir^ candidates, 
with D — > Kg7r + 7r~ , divided between the (a,b) long and (c,d) downstream K® categories. Fit 
results, including the signal and background components, are superimposed. 
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Figure 3: Invariant mass distributions of (a) B^ — > DK^ and (b) B^ — > Dn^ candidates, with 
D -)• K®K + K~ , shown with both K° cate gories combined. Fit results, including the signal and 
background components, are superimposed. 



used in a similar analysis of B ± — > DK ± decays, with D — > K ±r K T , K + K and tt + tt [T]. 
In that analysis the shape was constructed by applying the selection to a large simulated 
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sample containing many common backgrounds, each weighted by its production rate and 
branching fraction. The invariant mass distribution for the surviving candidates was cor- 
rected to account for small differences in resolution and PID performance between data 
and simulation, and two background PDFs were extracted by kernel estimation 25 ; one 
for B ± —7- DK ± and one for B — > Dn^ 1 decays. The partially reconstructed background 
PDFs are found to give a good description of both K® categories. 

An additional and significant background component exists in the B — > DK ± sample, 
arising from the dominant B ± — > Dn^ 1 decay on those occasions where the bachelor 
particle is misidentified as a kaon by the RICH system. In contrast, the B ± — > DK ± 
contamination in the B ± — > sample can be neglected. The size of this background is 
calculated through knowledge of PID and misidentification efficiencies, which are obtained 
from large samples of kinematically selected D*^ — > Dir^, D — > K^-r^ decays. The 
kinematic properties of the particles in the calibration sample are reweighted to match 
those of the bachelor particles in the B decay sample, thereby ensuring that the measured 
PID performance is representative of that in the B decay sample. The efficiency to 
identify a kaon correctly is found to be around 86%, and that for a pion to be around 
96%. The misidentification efficiencies are the complements of these numbers. From this 
information and from knowledge of the number of reconstructed B — > Dn^ 1 decays, the 
amount of this background surviving the B ± — > DK^ selection can be determined. The 
invariant mass distribution of the misidentified candidates is described by a Crystal Ball 



function 26 with the tail on the high mass side, the parameters of which are fitted in 
common between all the B ± — > DK ± samples. 

The number of B^ — > DK ± candidates in all categories is determined by TZ, and the 
number of — > Dn^ 1 events in the corresponding category. The ratio 1Z is determined 
in the fit and measured to be 0.085±0.005 (statistical uncertainty only) and is consistent 
with that observed in Ref. [I]. The yields returned by the invariant mass fit in the full 
fit region are scaled to the signal region, defined as 5247-5317 MeV/c 2 , and are presented 
in Tables [I] and [2] for the D — > K® tt + it~ and D — > K®K + K~ selections respectively. In 
the B — > D(KgTi + 'K~)K ± sample there are 654 ± 28 signal candidates, with a purity of 
86%. The corresponding numbers for the B ± — > D(K®K + K~)K ± sample are 102 ±5 and 
88%, respectively. The contamination in the B ± — > DK ± selection receives approximately 
equal contributions from misidentified B ± — > decays, combinatoric background and 
partially reconstructed decays. The partially reconstructed component in the signal region 
is dominated by decays of the type B — > Dp, in which a charged pion from the p decay is 
misidentified as the bachelor kaon, and B — > D*^, again with a misidentified pion. 

The Dalitz plots for B — > DK ± data in the signal region for the two D ->• K^h+h- 
final states are shown in Fig. |4j Separate plots are shown for B + and B~ decays. 

5 Binned Dalitz fit 

The purpose of the binned Dalitz plot fit is to measure the CP-violating parameters x± 
and y±, as introduced in Sect. [2j Following Eq. ^ these parameters can be determined 
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Table 1: Yields and statistical uncertainties in the signal region from the invariant mass fit, 
scaled from the full fit mass range, for candidates passing the B ± -»• Dh ± , D -> K s °7r+7r- 
selection. Values are shown separately for candidates containing long and downstream K® 
decays. The signal region is between 5247 MeV/c 2 and 5317 MeV/c 2 and the full fit range is 
between 5110 MeV/c 2 and 5800 MeV/c 2 . 





_B ± ->■ DK*- selection 


B ± ->■ Dtt 


± selection 


Fit component 


Lonj 


y 


Downstream 


Long 


Downstream 


B ± ->■ DK* 


213 ± 


13 


441 ± 25 






B ± ->■ Dit± 


11 ± 


3 


22 ±5 


2809 ± 56 


5755 ± 82 


Combinatoric 


9± 


4 


29 ±6 


22 ±3 


90 ±7 


Partially reconstructed 


11 ± 


1 


25 ±2 


25 ±1 


55 ± 1 



Table 2: Yields and statistical uncertainties in the signal region from the invariant mass fit, 
scaled from the full fit mass range, for candidates passing the — > Dh^, D — > K®K + K~ 
selection. Values are shown separately for candidates containing long and downstream K® 
decays. The signal region is between 5247 MeV/c 2 and 5317 MeV/c 2 and the full fit range is 
between 5110 MeV/c 2 and 5800 MeV/c 2 . 





B ± ->■ selection 


B ± — > selection 


Fit component 


Long 


Downstream 


Long 


Downstream 


B ± ->■ DK* 


32 ±2 


70 ±4 






B ± ->• Dtt ± 


1.6 ± 1.2 


3.4 ± 1.8 


417 ±20 


913 ±29 


Combinatoric 


0.6 ±0.5 


2.5 ±0.9 


4.8 ± 1.4 


18 ±2 


Partially reconstructed 


2.2 ±0.4 


2.9 ±0.5 


3.7 ±0.3 


7.7 ±0.5 



from the populations of each B ± — > DK ± Dalitz plot bin given the external information 
that is available for the q, and Ki parameters. In order to know the signal population 
in each bin it is necessary both to subtract background and to correct for acceptance 
losses from the trigger, reconstruction and selection. 

Although the absolute numbers of B + and B~ decays integrated over the Dalitz plot 
have some dependence on x± and y±, the additional sensitivity gained compared to using 
just the relative bin-to-bin yields is negligible, and is therefore not used. Consequently 
the analysis is insensitive to any B production asymmetries, and only knowledge of the 
relative acceptance is required. The relative acceptance is determined from the control 
channel B ± — > Dtt^. In this decay the ratio of b — > ucd to b — > cud amplitudes is expected 
to be very small (~ 0.005) and thus, to a good approximation, interference between the 
transitions can be neglected. Hence the relative population of decays expected in each 
B ± ->■ Dtt ± Dalitz plot bin can be predicted using the Ki values calculated with the 
D — > K®h + h~ model. Dividing the background-subtracted yield observed in each bin by 
this prediction enables the relative acceptance to be determined, and then applied to the 
B ± — > DK ± data. In order to optimise the statistical precision of this procedure, the bins 
+i and —i are combined in the calculation, since the efficiencies in these symmetric regions 
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Figure 4: Dalitz plots of — > DK^ candidates in the signal region for (a,b) D — > K®ir + -K~ 
and (c,d) D — > K®K + K~ decays, divided between (a,c) B + and (b,d) B~ . The boundaries of 
the kinematically-allowed regions are also shown. 



are expected to be the same in the limit that there are no charge-dependent reconstruction 
asymmetries. It is found that the variation in relative acceptance between non-symmetric 
bins is at most ~ 50%, with the lowest efficiency occurring in those regions where one of 
the pions has low momentum. 

Separate fits are performed to the B + and B~ data. Each fit simultaneously considers 
the two categories, the B ± — > DK ± and B ± — > Dn^ 1 candidates, and the two D — > 
K®h + h~ final states. In order to assess the impact of the D — > K®K + K~ data the fit is 
then repeated including only the D — > i\^7r + 7r~ sample. The PDF parameters for both 
the signal and background invariant mass distributions are fixed to the values determined 
in the global fit. The yields of all the background contributions in each bin are free 
parameters, apart from bins where a very low contribution is determined from an initial 
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Table 3: Results for x± and y± from the fits to the data in the case when both D — > 
Kg7r + 7r~ and D — > K®K + K~ are considered and when only the D — > Kg-K + 7i^ final state 
is included. The first, second, and third uncertainties are the statistical, the experimental 
systematic, and the error associated with the precision of the strong-phase parameters, 
respectively. The correlation coefficients are calculated including all sources of uncertainty 
(the values in parentheses correspond to the case where only the statistical uncertainties 
are considered). 



Parameter 


All data 


D i^s7r + 7r- alone 




xlO- 2 ] 


0.0 ± 4.3 ± 1.5 ±0.6 


1.6 ±4.8 ±1.4 ±0.8 


y- 


xlO- 2 ] 


2.7 ±5.2 ±0.8 ±2.3 


1.4 ±5.4 ±0.8 ±2.4 


corr 


[x-,y-) 


-0.10 (-0.11) 


-0.12 (-0.12) 


x + 


xlO- 2 ] 


-10.3 ±4.5 ±1.8 ±1.4 


-8.6 ±5.4 ±1.7 ±1.6 


y + 


xlO" 2 ] 


-0.9 ±3.7 ±0.8 ±3.0 


-0.3 ±3.7 ±0.9 ±2.7 


corr 




0.22 (0.17) 


0.20 (0.17) 



fit, in which case they are fixed to zero, to facilitate the calculation of the error matrix. 
The yields of signal candidates for each bin in the B ± — > Dn^ 1 sample are also free 
parameters. The amount of signal in each bin for the B — > DK ± sample is determined 
by varying the integrated yield and the x± and y± parameters. 

A large ensemble of simulated experiments are performed to validate the fit procedure. 
In each experiment the number and distribution of signal and background candidates are 
generated according to the expected distribution in data, and the full fit procedure is 
then executed. The values for x± and y± are set close to those determined by previous 



measurements 14 . It is found from this exercise that the errors are well estimated. 
Small biases are, however, observed in the central values returned by the fit and these 
are applied as corrections to the results obtained on data. The bias is (0.2 — 0.3) x 10~ 2 
for most parameters but rises to 1.0 x 10~ 2 for y + . This bias is due to the low yields in 
some of the bins and is an inherent feature of the maximum likelihood fit. This behaviour 
is associated with the size of data set being fit, since when simulated experiments are 
performed with larger sample sizes the biases are observed to reduce. 

The results of the fits are presented in Table [3] The systematic uncertainties are 
discussed in Sect. [6j The statistical uncertainties are compatible with those predicted by 
simulated experiments. The inclusion of the D — > K®K + K~ data improves the precision 
on x± by around 10%, and has little impact on y±. This behaviour is expected, as the 
measured values of Cj in this mode, which multiply x± in Eq. Q, are significantly larger 
than those of Sj, which multiply y±. The two sets of results are compatible within the 
statistical and uncorrelated systematic uncertainties. 

The measured values of (x±, y±) from the fit to all data, with their statistical likelihood 
contours are shown in Fig. |5j The expected signature for a sample that exhibits CP- 
violation is that the two vectors defined by the coordinates (x_,y_) and (x + ,y + ) should 
both be non-zero in magnitude, and have different phases. The data show this behaviour, 
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Figure 5: One (solid), two (dashed) and three (dotted) standard deviation confidence levels for 
(x+,y+) (blue) and (x_,y_) (red) as measured in — > DK^ decays (statistical only). The 
points represent the best fit central values. 



but are also compatible with the no CP violation hypothesis. 

In order to investigate whether the binned fit gives an adequate description of the 
data, a study is performed to compare the observed number of signal candidates in each 
bin with that expected given the fitted total yield and values of x± and y±. The number of 
signal candidates is determined by fitting in each bin for the B ± — > DK ± contribution for 
long and downstream K® decays combined, with no assumption on how this component 
is distributed over the Dalitz plot. Figure [6] shows the results in effective bin number 
separately for Nb++b~, the sum of B + and B~ candidates, which is a CP-conserving 
observable, and for the difference Nb+-b- > which is sensitive to CP violation. The effective 
bin number is equal to the normal bin number for B + , but is defined to be this number 
multiplied by —1 for B~ . The expectations from the (x±, y±) fit are superimposed as is, 
for the N b +-b- distribution, the prediction for the case x± = y± = 0. Note that the zero 
CP violation prediction is not a horizontal line at N b +~b- — because it is calculated 
using the total B + and B~ yields from the full fit, and using bin efficiencies that are 
determined separately for each sample. The data and fit expectations are compatible for 
both distributions yielding a \ 2 probability of 10% for N b ++b~ and 34% for N b +-b-- 
The results for the Nb+-b- distribution are also compatible with the no CP- violation 
hypothesis (x 2 probability = 16% ). 
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Figure 6: Signal yield in effective bins compared with prediction of (x±,y±) fit (black histogram) 
for D — > K®ir + 7T~ and D — > K®K + K~. Figure (a) shows the sum of B + and B~ yields. Figure 
(b) shows the difference of B + and B~ yields. Also shown (dashed line and grey shading) is the 
expectation and uncertainty for the zero CP-violation hypothesis. 



6 Systematic uncertainties 

Systematic uncertainties are evaluated for the fits to the full data sample and are presented 
in Table |4j In order to understand the impact of the CLEO-c (cj,Sj) measurements the 
errors arising from this source are kept separate from the other experimental uncertainties. 
Table[5]shows the uncertainties for the case where only D — > fC?7r + 7r~ decays are included. 
Each contribution to the systematic uncertainties is now discussed in turn. 

The uncertainties on the shape parameters of the invariant mass distributions as de- 
termined from the global fit when propagated through to the binned analysis induce 
uncertainties on x± and y±. In addition, consideration is given to certain assumptions 
made in the fit. For example, the slope of the combinatoric background in the data set 
containing D — > K®K + K~ decays is fixed to be zero on account of the limited sample 
size. The induced errors associated with these assumptions are evaluated and found to 
be small compared to those coming from the parameter uncertainties themselves, which 
vary between 0.4 x 10 -2 and 0.6 x 10~ 2 for the fit to the full data sample. 

The analysis assumes an efficiency that is flat across each Dalitz plot bin. In reality 
the efficiency varies, and this leads to a potential bias in the determination of x± and 
y±, since the non-uniform acceptance means that the values of (q,Sj) appropriate for 
the analysis can differ from those corresponding to the flat-efficiency case. The possible 
size of this effect is evaluated in LHCb simulation by dividing each Dalitz plot bin into 
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Table 4: Summary of statistical, experimental and strong-phase uncertainties on x± and y± 
in the case where both D — > K®tt + tt~ - and D K®K + K~ decays are included in the fit. All 
entries are given in multiples of 10~ 2 . 



Component 


0"(x_) 




a(x+) 




Statistical 


4.3 


5.2 


4.5 


3.7 


VJlUUdl lit ollcLUc UcLl cLlllc L"l o 


n 4 


n 4 


u.u 


n 4 


Efficiency effects 


0.3 


0.4 


0.3 


0.4 


CP violation in control mode 


1.3 


0.4 


1.5 


0.2 


Migration 


0.4 


0.2 


0.4 


0.2 


Partially reconstructed background 


0.2 


0.3 


0.2 


0.2 


PID efficiency 


0.1 


0.2 


0.2 


< 0.1 


Shape of misidentified 5 ± — > Dir^ 1 


0.1 


0.1 


0.3 


< 0.1 


Bias correction 


0.2 


0.3 


0.2 


0.5 


Total experimental systematic 


1.5 


0.9 


1.8 


0.8 


Strong-phase systematic 


0.6 


2.3 


1.4 


3.0 



Table 5: Summary of statistical, experimental and strong-phase uncertainties on x± and y± 
in the case where only D — > Kg-K + 7r~ decays are included in the fit. All entries are given in 
multiples of 10 -2 . 



Component 


<j(x_) 


o-(y-) 






Statistical 


4.8 


5.4 


5.4 


3.7 


Global fit shape parameters 


0.4 


0.4 


0.6 


0.4 


Efficiency effects 


0.2 


0.2 


0.3 


0.4 


CP violation in control mode 


1.2 


0.5 


1.5 


0.2 


Migration 


0.4 


0.2 


0.4 


0.2 


Partially reconstructed background 


0.1 


0.1 


0.3 


0.2 


PID efficiency 


< 0.1 


0.2 


< 0.1 


< 0.1 


Shape of misidentified — > Dir^ 1 


0.1 


< 0.1 


0.1 


< 0.1 


Bias correction 


0.2 


0.3 


0.2 


0.6 


Total experimental systematic 


1.4 


0.8 


1.7 


0.9 


Strong-phase systematic 


0.8 


2.4 


1.6 


2.7 



many smaller cells, and using the BaBar amplitude model [5j[6] to calculate the values 
of Ci and Sj within each cell. These values are then averaged together, weighted by the 
population of each cell after efficiency losses, to obtain an effective (cj, Si) for the bin 
as a whole, and the results compared with those determined assuming a flat efficiency. 
The differences between the two sets of results are found to be small compared with the 
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CLEO-c measurement uncertainties. The data fit is then rerun many times, and the 
input values of (cj, s») are smeared according to the size of these differences, and the mean 
shifts are assigned as a systematic uncertainty. These shifts vary between 0.2 x 10 -2 and 
0.3 x 10" 2 . 

The relative efficiency in each Dalitz plot bin is determined from the B — > 
control sample. Biases can enter the measurement if there are differences in the relative 
acceptance over the Dalitz plot between the control sample and that of signal B — > DK ± 
decays. Simulation studies show that the acceptance shapes are very similar between the 
two decays, but small variations exist which can be attributed to kinematic correlations 
induced by the different PID requirements on the bachelor particle from the B decay. 
When included in the data fit, these variations induce biases that vary between 0.1 x 10~ 2 
and 0.3 x 10~ 2 . In addition, a check is performed in which the control sample is fitted 
without combining together bins +i and — i in the efficiency calculation. As a result of this 
study small uncertainties of < 0.3 x 10~ 2 are assigned for the D — > K®K + K~ measurement 
to account for possible biases induced by the difference in interaction cross-section for K~ 
and K + mesons interacting with the detector material. These contributions are combined 
together with the uncertainty arising from efficiency variation within a Dalitz plot bin to 
give the component labelled 'Efficiency effects' in Tables [4] and [5j 

The use of the control channel to determine the relative efficiency on the Dalitz plot 
assumes that the amplitude of the suppressed tree diagram is negligible. If this is not the 
case then the B~ final state will receive a contribution from D° decays, and this will lead 
to the presence of CP violation via the same mechanism as in B — > DK decays. The size 
of any CP violation that exists in this channel is governed by r^ n , 7 and 5^, where the 
parameters with superscripts are analogous to their counterparts in B — > DK ± decays. 
The naive expectation is that 0.005 but larger values are possible, and the studies 

reported in Ref. [l] are compatible with this possibility. Therefore simulated experiments 
are performed with finite CP violation injected in the control channel, conservatively set- 
ting r^ 71 " to be 0.02, taking a wide variation in the value of the unknown strong-phase 
difference 5^, and choosing 7 = 70°. The experiments are fit under the no CP violation 
hypothesis and the largest shifts observed are assigned as a systematic uncertainty. This 
contribution is the largest source of experimental systematic uncertainty in the measure- 
ment, for example contributing an error of 1.5 x 10~ 2 in the case of x + in the full data 
fit. 

The resolution of each decay on the Dalitz plot is approximately 0.004 GeV 2 /c 4 for 
candidates with long K® decays and 0.006 GeV 2 /c 4 for those containing downstream K® 
in the m 2 + and m 2 directions. This is small compared to the typical width of a bin, 
nonetheless some net migration is possible away from the more densely populated bins. 
At first order this effect is accounted for by use of the control channel, but residual effects 
enter because of the different distribution in the Dalitz plot of the signal events. Once 
more a series of simulated experiments is performed to assess the size of any possible bias 
which is found to vary between 0.2 x 10" 2 and 0.4 x 10" 2 . 

The distribution of the partially reconstructed background is varied over the Dalitz 
plot according to the uncertainty in the make-up of this background component. From 
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these studies an uncertainty of (0.2 — 0.3) x 10~ 2 is assigned to the fit parameters in the 
full data fit. 

Two systematic uncertainties are evaluated that are associated with the misidentified 
B — > Dtt^ 1 background in the B — > DK ± sample. Firstly, there is a 0.2 x 10~ 2 
uncertainty on the knowledge of the efficiency of the PID cut that distinguishes pions 
from kaons. This is found to have only a small effect on the measured values of x± 
and y±. Secondly, it is possible that the invariant mass distribution of the misidentified 
background is not constant over the Dalitz plot, as is assumed in the fit. This can occur 
through kinematic correlations between the reconstruction efficiency on the Dalitz plot 
of the D decay and the momentum of the bachelor pion from the B ± decay. Simulated 
experiments are performed with different shapes input according to the Dalitz plot bin 
and the results of simulation studies, and these experiments are then fitted assuming a 
uniform shape, as in data. Uncertainties are assigned in the range (0.1 — 0.3) x 10~ 2 . 

An uncertainty is assigned to each parameter to accompany the correction that is 
applied for the small bias which is present in the fit procedure. These uncertainties are 
determined by performing sets of simulated experiments, in each of which different values 
of x± and y± are input, corresponding to a range that is wide compared to the current 
experimental knowledge, and also encompassing the results of this analysis. The spread 
in observed bias is taken as the systematic error, and is largest for y + , reaching a value 
of 0.5 x 10~ 2 in the full data fit. 

Finally, several robustness checks are conducted to assess the stability of the re- 
sults. These include repeating the analysis with alternative binning schemes for the 
D — > A'g7r + 7r~ data and performing the fits without making any distinction between K® 
category. These tests return results compatible with the baseline procedure. 

The total experimental systematic uncertainty from LHCb-related sources is deter- 
mined to be 1.5 x 10~ 2 on x_, 0.9 x 10 -2 on 1.8 x 10~ 2 on x + and 0.8 x 10~ 2 on 
y + . These are all smaller than the corresponding statistical uncertainties. The dominant 
contribution arises from allowing for the possibility of CP violation in the control channel, 
B —7- Dir. In the future, when larger data sets are analysed, alternative analysis methods 
will be explored to eliminate this potential source of bias. 

The limited precision on (q, Sj) coming from the CLEO-c measurement induces un- 
certainties on x± and y± [12]. These uncertainties are evaluated by rerunning the data 
fit many times, and smearing the input values of (cj, Sj) according to their measurement 
errors and correlations. Values of (0.6 — 3.0) x 10 -2 are found for the fit to the full sample. 
When evaluated for the D — > K® r n +r n~ data set alone, the results are similar in magnitude, 



but not identical, to those reported in the corresponding Belle analysis 13 . Differences 
are to be expected, as these uncertainties have a dependence on the central values of the 
x± and y± parameters, and are sample-dependent for small data sets. Simulation studies 
indicate that these uncertainties will be reduced when larger B — > DK ± data sets are 
analysed. 

After taking account of all sources of uncertainty the correlation coefficient between 
X- and y_ in the full fit is calculated to be —0.10 and that between x + and y + to be 
0.22. The correlations between B~ and B + parameters are found to be small and can 
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be neglected. These correlations are summarised in Table [3j together with those coming 
from the statistical uncertainties alone, and those from the fit to D — > K®tt + it~ data. 



7 Interpretation 



The results for x± and y± can be interpreted in terms of the underlying physics param- 
eters 7, tb and 5b- This is done using a frequentist approach with Feldman-Cousins 
ordering [27], using the same procedure as described in Ref. (l3l. In this manner confi- 



using the same procedure as described in Ref. 
dence levels are obtained for the three physics parameters. The confidence levels for one, 
two and three standard deviations are taken at 20%, 74% and 97%, which is appropriate 
for a three-dimensional Gaussian distribution. The projections of the three-dimensional 
surfaces bounding the one, two and three standard deviation volumes onto the (7, tb) 
and (7, 5b) planes are shown in Fig. [7J The LHCb-related systematic uncertainties are 
taken as uncorrelated and correlations of the CLEO-c and statistical uncertainties are 
taken into account. The statistical and systematic uncertainties on x and y are combined 
in quadrature. 
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Figure 7: Two-dimensional projections of confidence regions onto the (7, re) and (7, #b) planes 
showing the one (solid) and two (dashed) standard deviations with all uncertainties included. 
For the (7, r#) projection the three (dotted) standard deviation contour is also shown. The 
points mark the central values. 



The solution for the physics parameters has a two-fold ambiguity, (7, 5b) and (7 + 
180°, 5 B + 180°). Choosing the solution that satisfies < 7 < 180° yields r B = 0.07±0.04, 
7 = (44lgg) and 5b = (137jj[g) . The value for r B is consistent with, but lower than, 
the world average of results from previous experiments 15 . This low value means that 



it is not possible to use the results of this analysis, in isolation, to set strong constraints 
on the values of 7 and 5b, as can be seen by the large uncertainties on these parameters. 



17 



8 Conclusions 

Approximately 800 B ± — > DK ± decay candidates, with the D meson decaying either 
to .K"g7r + 7r~ or K®K + K~ , have been selected from 1.0 fb _1 of data collected by LHCb 
in 2011. These samples have been analysed to determine the CP-violating parameters 
x± = r B cos(5b ± 7) and y± = r B sin(5 B ± 7), where r B is the ratio of the absolute values 
of the B + — > D°K~ and B + — > D°K~ amplitudes, 5 B is the strong-phase difference 
between them, and 7 is the angle of the unitarity triangle. The analysis is performed 
in bins of D decay Dalitz space and existing measurements of the CLEO-c experiment 
are used to provide input on the D decay strong-phase parameters (cj,Sj) [12]. Such 
an approach allows the analysis to be essentially independent of any model-dependent 
assumptions on the strong phase variation across Dalitz space. It is the first time this 
method has been applied to D — > K®K + K~ decays. The following results are obtained 

= (0.0 ± 4.3 ± 1.5 ± 0.6) x 10" 2 , y_ = (2.7 ± 5.2 ± 0.8 ± 2.3) x 10~ 2 , 

x+ = (-10.3 ± 4.5 ± 1.8 ± 1.4) x 10~ 2 , y + = (-0.9 ± 3.7 ± 0.8 ± 3.0) x 10~ 2 , 



where the first uncertainty is statistical, the second is systematic and the third arises from 
the experimental knowledge of the (cj, Sj) parameters. These values have similar precision 
to those obtained in a recent binned study by the Belle experiment [l3| . 

When interpreting these results in terms of the underlying physics parameters it is 
found that r B = 0.07 ± 0.04, 7 = (44±f|) and 5 B = (137^)°. These values are 



consistent with the world average of results from previous measurements 15 , although the 



uncertainties on 7 and S B are large. This is partly driven by the relatively low central value 
that is obtained for the parameter r B . More stringent constraints are expected when these 
results are combined with other measurements from LHCb which have complementary 
sensitivity to the same physics parameters. 
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